Elliptic Flow at SPS and RHIC: From Kinetic Transport to Hydrodynamics 



P.F. Kolb Q - c , P. Huovinen d , U. Heinz a - C , and H. Heiselberg 6 
a Theoretical Physics Division, CERN, CH-1211 Geneva 23, Switzerland 
b Institut fur Theoretische Physik, Universitat Regensburg, D-93040 Regensburg, Germany 
c Department of Physics, The Ohio State University, 174 West 18th Ave., Columbus, OH 43210, USA 
d Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 
e NORDITA, Blegdamsvej 17, DK-2100 Copenhagen 0, Denmark 
(February 1, 2008) 



O 
O 

o 

O 

Q 



(N 
> 

m 

(N 

o 
o 

Oh! 

i 

^ ■ 



X 



Anisotropic transverse flow is studied in Pb+Pb and 
Au+Au collisions at SPS and RHIC energies. The centrality 
and transverse momentum dependence at midrapidity of the 
elliptic flow coefficient i>2 is calculated in the hydrodynamic 
and low density limits. Hydrodynamics is found to agree well 
with the RHIC data for semicentral collisions up to transverse 
momenta of 1—1.5 GeV/c, but it considerably overestimates 
the measured elliptic flow at SPS energies. The low density 
limit LDL is inconsistent with the measured magnitude of i>2 
at RHIC energies and with the shape of its pt-dependence at 
both RHIC and SPS energies. The success of the hydrody- 
namic model points to very rapid thermalization in Au+Au 
collisions at RHIC and provides a serious challenge for kinetic 
approaches based on classical scattering of on-shell particles. 

PACS numbers: 25.75-q, 24.85. +p, 25.75.Ld 
Keywords: Relativistic heavy-ion collisions; elliptic 
flow 

1. Introduction.- Anisotropic flow has been measured 
in relativistic nuclear collisions and calculations, in 

particular for collisions at SPS and RHIC energies, exist 
within a variety of frameworks: hydrodynamics [p|-^0[, 
the low density limit of kinetic theory p arton cas- 
cades Jl2| , |l3[ , hadronic cascade cod es |l4] ~ 19[ , combina- 
tions thereof |2(J, and jet quenching ||2lf~Recent detailed 
SPS HI] and RHIC H data on the centrality and trans- 
verse momentum dependence of pion and proton elliptic 
flow begin to discriminate between different model cal- 
culations and to yield quantitative insights into initial 
conditions, compression, rescattering time scales, expan- 
sion dynamics and possible phase transitions during the 
expansion stage of the reaction zone. 

Cascade calculations based on the incoherent scatter- 
ing of classical on-shell particles and the low density limit 
of classical kinetic theory are expected to work best for 
peripheral collisions where the density of produced parti- 
cles is sufficiently low and only a few rescatterings occur. 
Central collisions produce higher particle densities where 
the hydrodynamic limit may be more suitable. One of the 
most interesting questions in the kinetic theory of rela- 
tivistic heavy ion collisions is where and how the transi- 
tion between the dilute and dense limits happens. 

In this paper we present detailed calculations of the 
impact parameter and transverse momentum dependence 
of elliptic flow in a hydrodynamic model [Qj8| and in the 
low density limit (LDL) of Ref. and we compare the 
results to recent SPS and RHIC data We test the 



validity of these two models in semi-central Pb+Pb and 
Au+Au collisions where the elliptic flow signal is large 
and thus presents a perfect tool. 

2. Geometry and Anisotropic Flow. - Consider two nu- 
clei of radius R colliding at impact parameter 6. We refer 
to the collision as "central" when b^^R, "semi-central" 
when and "peripheral" when In 

non-central collisions the initial overlap zone is narrower 
parallel than perpendicular to the impact parameter. A 
simple ansatz for the initial density in the interaction 
region is that it scales with the number of participat- 
ing nucleons per unit area in the transverse plane 
This predicts a nearly linear dependence of the multiplic- 
ity on ./Vpart in central and semi-central Au+Au collisions 
at RHIC p2fl . Other initializations have been proposed 
p3[ . The sensitivity of elliptic flow to the details of the 
initialization is studied elsewhere |22| . 




FIG. 1. The number of participants Ap art , the initial trans- 
verse radii R x , R y , and the initial spatial deformation e x = S 
for Au+Au at ^/s = 130^4GeV, as functions of the scaled 
impact parameter b/R. For details see text. 

The nuclear thickness function entering the calculation 
of the initial conditions is evaluated with realistic 

Woods-Saxon nuclear density profiles for the incoming 
nuclei, with radius R = 1.12 A x / 3 -0.86 A" 1 / 3 fm (A=207 
for Pb+Pb at the SPS and A=197 for Au+Au at RHIC) 
and standard surface diffuseness a — 0.54 fm Jt|]. The 
nucleon-nucleon cross section is taken as 32 mb at y/s — 
17GeV and 40 mb at V? = 130 GeV. 

With the same parametrization of the initial trans- 
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verse density profile we compute the initial transverse 
source radii parallel (R x (b) = (a; 2 ) 1 / 2 ) and perpendicu- 
lar (R y (b) = (y 2 ) 1 / 2 ) to the impact parameter. These as 
well as the initial spatial deformation 



density has no measurable consequences for the develop- 
ing flow pattern, since the pressure is insensitive to nj, 
when the latter is small M. 
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are shown in Fig. |l] for Au+Au at yfs = 130 A GeV as 
functions of the scaled impact parameter b/R. 

Particles are initially produced azimuthally symmetric 
in momentum space. The initial deformation in coor- 
dinate space generates an azimuthally asymmetric mo- 
mentum distribution if and only if the produced parti- 
cles rescatter off each other. The more rescatterings take 
place and the larger the initial deformation <5, the more 
anisotropic the final momentum distribution will be. A 
quantitative measure of this anisotropy is provided by its 
harmonic coefficients v n in an event-by-event Fourier ex- 
pansion with respect to the azimuthal angle </> |24| . v-i is 
called the "elliptic flow coefficient" . 

To relate the initial spatial anisotropy 5 to the mea- 
sured momentum anisotropy vi one has to model the 
interactions during the fireball expansion. The initial 
average density and source size is larger in central than 
in peripheral collisions. The particle mean free paths 
A m fp = 1/ ipp) i where a is the scattering cross section and 
p the time-dependent average density of scatterers, vary 
from a few fm or less initially to infinity at freeze-out. 
Since the ratio \ m { p /R Xt y is smaller over a longer time in 
central than in peripheral collisions, hydrodynamic mod- 
els are more likely to work in central or semicentral colli- 
sions, especially at high energies where the initial particle 
densities are large. The low density limit (LDL), on the 
other hand, should apply to peripheral collisions, espe- 
cially at low energies with small initial particle densities. 

3. The Hydrodynamic Limit. - The full hydrodynamic 
treatment of a non-central collision is a tedious 3+1 di- 
mensional problem |^5],[l0 26 1. We reduce the complex- 
ity of the task to 2+1 dimensions by assuming boost- 
invariant longitudinal flow ^ 0]. This assumption lim- 
its our description to a region around midrapidity which 
is, however, expected to grow as the collision energy in- 
creases. 

The evolution of a hydrodynamical system is deter- 
mined by its initial conditions and equation of state 
(EOS). We fix the initial conditions as in jjJD by requir- 
ing a good fit to the p t -spectra of protons and negatively 
charged particles in central Pb+Pb collisions at the SPS 
0. The SPS initial conditions are scaled to the RHIC 
energy of s/s = 130 A GeV by adjusting the initial energy 
density eo (keeping the product Tqtq of the initial tem- 
perature and thermalization time fixed Q|) until the final 
charged particle pseudorapidity density at midrapidity 
agrees with the published measurement by the PHOBOS 
Collaboration [£7j . The initial baryon density was chosen 
to give the ratio p/p = 0.65 in the final state pS| ]. Ad- 
justing the initial baryon density at fixed initial energy 
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Table 1. Freeze-out temperatures, equations of state and 
initial conditions for central (b = 0) collisions employed 
for the hydrodynamical calculations shown in this paper. 
The last two rows show the final charged particle mul- 
tiplicity densities in rapidity and pseudorapidity for the 
6% most central collisions, to facilitate comparison with 
the PHOBOS data ||. 

We use two different equations of state to check how 
the quark-hadron phase transition or its absence affects 
the flow anisotropy. EOS Q has a first order phase tran- 
sition to QGP at T c = 165 MeV whereas EOS H contains 
only hadronic resonances at all densities. At RHIC ener- 
gies EOS H leads to unrealistic particle densities in the 
initial state. However, this EOS gives an impression how 
a system without a phase transition would behave; this 
is important when trying to separate specific phase tran- 
sition signatures |?],|8| from generic hydrodynamical fea- 
tures. Further details of the construction of these equa- 
tions of state can be found in (2{| . 

There is still no final consensus at which temperature 
the hydrodynamical description breaks down at SPS en- 
ergies pot , nor is there any a priori reason why this tem- 
perature should be the same at SPS and RHIC. In our 
earlier studies of Pb+Pb collisions at the SPS Q we 
found EOS Q to favour a lower freeze-out temperature 
(T f w 120 MeV) than EOS H (T f w 140 MeV). Still, a 
fit of the pion and proton p t spectra from 158 A GeV 
Pb+Pb collisions H with EOS H and T f w 120 MeV 
looks quite acceptable |33[ . Therefore we do the SPS 
calculations using EOS Q with freeze-out energy density 
e £ = 0.06 GeV/fm 3 (corresponding to T f « 120 MeV), 
and we probe the effect of the freeze-out temperature 
on elliptic flow by using EOS H with freeze-out energy 
densities ef = 0.06 as well as 0.15 GeV/fm 3 (the latter 
corresponds to Tf an 140 MeV). 

At RHIC the baryon number density at freeze-out is 
much smaller than at the SPS, and the energy densities 
corresponding to T { = 120 and 140 MeV are e f R* 0.05 
and 0.14 GeV/fm 3 , respectively. We do the RHIC cal- 
culations for EOS Q with both of these freeze-out en- 
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ergy densities and for EOS H with freeze-out density 
£f = 0.14 GeV/fm 3 . To convert the fluid variables to par- 
ticle and resonance distributions we employ the Cooper- 
Frye freeze-out prescription Subsequent resonance 

decays are calculated using the decay kinematics de- 
scribed in |p5| . The elliptic flow parameters i>2 (p*) and 
V2 are obtained by Fourier expanding the calculated dif- 
ferential and p t -integrated momentum distributions. For 
Tf pa 140 MeV decay contributions reduce the elliptic flow 
of pions by 18-25% depending on the EOS and impact 
parameter; the corresponding decrease at Tf ss 120 MeV 
is 8-15%f| 

The pt-averaged elliptic flow V2 as a function of col- 
lision centrality is shown in Fig. ^ for Pb+Pb at the 
SPS and in Fig. | for Au+Au at RHIC. At the SPS en- 
ergy the calculation is for midrapidity pions of all pt 
whereas the RHIC results include all charged particles 
with |r?| < 1.3 and 0.1 < p t < 2GeV/c f§. We have used 
the same centrality measures (impact parameter b at the 
SPS, the fractional charged particle multiplicity density 
at midrapidity Noh/N max at RHIC) as in the correspond- 
ing experimental publications.^] 

At the SPS energy the hydrodynamical calculation can 
reproduce the data only for the most central collisions. 
Already in semi-central collisions the calculations over- 
predict the measured elliptic flow significantly, and the 
disagreement increases to about a factor 2 in peripheral 
collisions. However, the comparison in Fig. || should be 
viewed with some care: as stated above, our calculation 
applies only to midrapidity (yi a b = 2.9) whereas the sta- 
tistically more significant open squares correspond to an 
average over the forward hemisphere. The preliminary 
data on the rapidity dependence of V2 reported in Q 
show non-trivial structures which are averaged out in the 
open squares. On the other hand, when taking into ac- 
count the p* cut in one finds that the midrapidity 
data from g at y = 3.25 (shown as crosses in Fig. g) 
are significantly lower than the extrapolation of v% to 
midrapidity published in ||. A reliable measurement of 
V2 at midrapidity is thus presently not available for this 



1 In |To| resonance decays were found to reduce the pion el- 
liptic flow by 40%. However, in that work only pions with 
50 < pt < 350 MeV/c were taken into account whereas we 
include either all p t (SPS) or only cut out p t < lOOMeV/c 
(RHIC). Resonance decays reduce the pion elliptic flow espe- 
cially at low values of p t ; using similar p t cuts as in ^| at 
SPS energies our results are compatible. 

2 To compare with the RHIC data, we had to scale iV m ax 
by a factor 0.95 since the hydrodynamic calculation cannot 
account for fluctuations in the charged multiplicity at fixed 
impact parameter: while our "wounded nucleon" initializa- 
tion nicely describes the measured multiplicity distribution 
M up to the "knee" at b = 0, the calculated distribution 
drops there abruptly to zero while the data show fluctuations 
up to a value -/V max which is ~ 5% larger. 



collision system, but would be very welcome. 
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FIG. 2. Elliptic flow for pions at midrapidity vs. centrality, 
for 158AGeV Pb+Pb collisions. Hydrodynamic calculations 
and results from the LDL are compared to NA49 data j3]Q. 
Numbers in brackets give Tf in MeV and the reduction from 
resonance decays, respectively. 

Taking the available data at face value, we are unable 
to account for them within the hydrodynamic approach 
even when varying the initial conditions, the EOS and the 
freeze-out temperature within the constraints provided 
by the measured single particle distributions. Several 
examples are shown in Fig. 0. 




FIG. 3. Centrality dependence of the elliptic flow coef- 
ficient «2 for charged particles from Au+Au collisions at 
y/~s = 130 A GeV. The data || are shown with the quoted 
systematic error of ±0.005. For details see text. 

Fig. U shows that, on the other hand, hydrodynam- 
ics successfully reproduces the elliptic flow measured at 
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RHIC for central and semi-central collisions. The best 
agreement is reached for EOS Q with Tf = 140MeV, 
where discrepancies begin to be significant only at im- 
pact parameters above 7 fm, but stay below 20% even for 
the most peripheral collisions. Lower freeze-out temper- 
atures or the use of EOS H (which is effectively harder 
at these collision energies, where the expanding matter 
spends a large fraction of its total lifetime in the soft 
phase transition region |^|) predict somewhat larger el- 
liptic flow. They also give larger radial flow; when single- 
particle spectra from this experiment become available, 
they will provide a crucial test of the approach and re- 
move the ambiguity of the freeze-out temperature |36|. 
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FIG. 4. The elliptic flow of charged particles from Au+Au 
collisions at ^/s = 130 A GeV vs. transverse momentum. Hy- 
drodynamic calculations and predictions from the LDL are 
compared with the data ||. The shape of the LDL curve 
reflects the weighting of different hadrons with their contri- 
bution to the charged particle spectrum; at small pt it is dom- 
inated by light pions, at high p t by heavy baryons. 



Fig. ^ shows the pt-dependence of elliptic flow for min- 
imum bias Au+Au collisions at RHIC [|| . Following the 
prescription in |5j it is calculated by 



Jbdb v 2 {p t ;b) 



dN cl 



dy p t dp, 



■(b) 



/•6d6 *v*_(6) 

J dy p t dp t v I 
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with a cut-off at i) max = 13.5 fm. The agreement between 
the data and the hydrodynamical calculations is excel- 
lent, especially when one considers the small variations 
of the latter upon changing parameters within the range 
allowed by the constraints. Only for p t above about 1.5 
GeV/c does the measured elliptic flow lag behind the 
hydrodynamic prediction, indicating a departure from 
thermalization for high-p t particles. In future it will be 
interesting to follow the data to higher p t where v 2 is 
first expected to saturate due to lack of thermalization, 
before decreasing again as expected from jet quenching 



pl| . The hydrodynamic curves start out quadratically at 
low p t , as required by general principles |37|] , then quickly 
p6| turn over to an approximately linear rise and keep 
increasing monotonically with p t , eventually saturating 
at v 2 (pt) = 1 as p t — > oo. 

It is not immediately obvious how the good agreement 
between theory and data in Fig. ^ is compatible with the 
clearly visible discrepancies in Fig. ^ for (semi-)peripheral 
collisions. One possibility is that the impact parame- 
ter dependence |22]| of the charged particle multiplicity 
(i.e. of the normalization of the spectra which enter the 
weighting procedure (|^)) is different in theory and ex- 
periment. This can be settled by measuring the depen- 
dence of dNch/dy on the number of participants iV par t 



||22| , [23|| 38] , and by providing v 2 (pt) for different central- 
ity bins. A more likely explanation is suggested by the 
observation that the clearly visible differences in Fig. || 
resulting from varying the EOS and Tf are much smaller 
in Fig. |j. An analysis shows that the variations in Fig. || 
result from different slopes of the single particle spectra 
which, when averaging the nearly identical hydrodynamic 
curves in Fig. ^ over p t , give different relative weights to 
the regions of small and large v 2 . The lower p t -averaged 
v 2 for peripheral collisions could thus be due to steeper 
single-particle spectra in peripheral collisions than pre- 
dicted by the model, e.g. due to earlier freeze-out at 
higher Tf and smaller radial flow (20) . This can be clari- 
fied by measuring the single-particle spectra [|6| and the 
Pt-dependence of elliptic flow at different centralities. 
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FIG. 5. pt-dependence of elliptic flow for pions (left) and 
protons (right) from Pb+Pb collisions at y/s = 17 A GeV with 
impact parameters b < 11 fm. The data Q correspond to 
6.5 fm < b <8fm and are averaged over the forward rapidity 
interval 4 < y < 5 while the hydrodynamic calculations apply 
to midrapidity y — 2.9. The squares show the NA49 data after 
correction for azimuthal HBT correlations J3!| . For details see 
text. 

The pt-dependence of v 2 for pions and protons from 
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semi-central Pb+Pb collisions at the SPS is shown in 
Fig. Unfortunately, no midrapidity data are available, 
and a meaningful comparison between theory and ex- 
periment of the magnitude of V2 is thus not possible. 
We note, however, that the data show the same approx- 
imately linear rise of vi with p t and the smaller elliptic 
flow at small pt for protons than for pions as predicted 
by hydrodynamics. Such a linear rise is inconsistent with 
the Low Density Limit (LDL) (see below); it remains to 
be seen to what extent realistic kinetic codes are able to 
reproduce a linear rise of i>2(Pt)> albeit with a smaller 
slope than predicted by hydrodynamics, at SPS energies. 

4- The Low Density Limit. - If the nuclear overlap zone 
is small and the initial density not very large, produced 
particles can escape from the reaction zone suffering only 
a few reinteractions. In the extreme limit the system is 
streaming freely, and no collective flow builds up at all. 
For sufficiently dilute systems the elliptic flow can be cal- 
culated from the first order correction to free streaming 
arising from particle collisions. Such a perturbative ap- 
proach is valid as long as the particle mean free paths 
A m f p are larger than the overlap zone R x , y (b). 

In the low density limit (LDL) the effect of scatterings 
was calculated to first order by inserting the free stream- 
ing distribution into the Boltzmann collision term JlT| . 
The resulting momentum space anisotropy leads to an 
elliptic flow coefficient v\ for particle species i. scattering 
with particles j (i,j — ir,p, K, . . .), given by [ [Ll| 



v\{b,pt) 



vL(Pt) 5(b) 



16tt R x {b)R y {b) 



dy 



Here, i>ij_(pt) is the transverse velocity of a particle i with 
momentum p t , and Vij is the relative velocity between 
particle i and scatterer j. Brackets (. . .) denote aver- 
aging over scatterer momenta pj. a^ T is the transport 
cross section responsible for momentum transfer. Sim- 
ilar to the hydrodynamic limit, the elliptic flow coeffi- 
cient is proportional to the initial spatial deformation. 
Its pt-dependence is, however, quite different, saturating 
at large pt (when Vix — * c) at values much below 1. At 
small p t , V2 rises quadratically as it should [jj7| , at a rate 
which is directly given by the asymptotic saturation level. 
This coupling between the curvature at small p t and the 
saturation level at large p t does not permit a sustained 
approximately linear rise of V2 in the intermediate pt- 
region. 

Let us now try to quantify the parameters entering the 
LDL formula (|) for SPS and RHIC energies. S(b) and 
R x ,y(b) are obtained from Fig. [|. We assume that in the 
region of applicability of the LDL the matter can be de- 
scribed with hadronic degrees of freedom. For isotropic 
scattering the transport cross sections are about half 
of the total cross section. However, the produced pi- 
ons are fast and scatter off other pions mostly in the 
forward direction; the same is true for p-wave p- and 
A-resonance scattering. We estimate that the angle- 
averaged cross sections are only ~ (sin 2 9 cos 9) = 1/4 



of the total ones, averaged over the relevant relative mo- 
mentum range. For these latter we assume cr 7rM ~ 10 nib, 
a 7rB = aP M _ 30 m b, and a pB ~ 40 mb where M and B 
stand for arbitrary mesonic and baryonic resonances. 

If the cross sections were the same for all particles 
we could replace the sum over dNj/dy by § x dN^/dy. 
We use dN^/dy — 410 in central Pb+Pb collisions at 
the SPS (Vs = 17,4GeV) @ and dN ch /dy ~ 630 for 
central Au+Au collisions at RHIC (y/s = 130^4GeV) 
p7p . [dN/dy at midrapidity is about 15% larger than 
dN/dr)l„ I < i .] Since the baryon and anti-baryon cross sec- 
tions are larger, their rapidity densities are required sepa- 
rately; we use dNj^ s /dy ~ 110 for central Pb+Pb colli- 
sions at the SPS Q and estimate roughly dN™"£ /dy ~ 
70 for central Au+Au collisions at RHIC. For the cen- 
trality dependence we assume that the rapidity densities 
scale with the number of participants as shown in Fig. [l]. 

All relative velocities are assumed to be of the order of 
the speed of light. To obtain the p t -averaged elliptic flow 
we average over with an exponential mt-distribution 
with inverse slope of 130 MeV; this gives (t)j_) = 0.68 



for pions and 



= 0.22 for protons. In the RHIC data 



only charged particles withp t > lOOMeV/c are included, 
which we try to take into account by increasing (v'j_) to 
0.75 (assuming pion dominance). 

Finally, we have to correct for resonance decay contri- 
butions. The actual rapidity density of scatterers in the 
reaction zone is smaller than the observed dN/dy since 
a large fraction of the latter comes from unstable reso- 
nances which only decay after the ^-generating rescat- 
terings have happened. We estimate the corresponding- 
reduction factor for the density of scatterers (and thus 
for V2) to be about a factor 2. For illustration of the 
systematic uncertainties we also show curves where only 
25% of the final charged multiplicity arise from resonance 
decays. Clearly, all these numbers are very rough, and 
our estimates for the factor multiplying vf ± in (ra) could 
be easily off by 50% in both directions. 

Curves showing the elliptic flow from the LDL are in- 
cluded in Figs. I |. The LDL gives about 50% more el- 
liptic flow at RHIC than at the SPS. This reflects mostly 
the corresponding ratio of the charged multiplicity den- 
sities as shown by the following expression (see Eq. (g|)): 



2 RHIC _ dN* mc /dy + crhic dN^f/dy 
dN^ s /dy + c SPS dN s B l s s /dy 



,,sps 



(4) 



The constant is given by c = (2/3)(er£ B /(7£ r M - 1) w 4/3. 
While the LDL thus happens to be able to reproduce 
the impact parameter dependence of the ^-averaged el- 
liptic flow at the SPS (Fig. ||), it slightly underpredicts 
the same quantity at RHIC (Fig. |J). Inspection of the 
pt-dependence of V2 in minimum bias events (which is 
dominated by semi-central collisions) reveals, however, 
that the LDL gets it completely wrong at RHIC energies 
(Fig. ||) . Fig. H shows that, while the pion flow data at the 
SPS are inconclusive due to their limited p t -coverage, the 
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Pt-dependence of v% for protons from semi-central Pb+Pb 
collisions at the SPS is again incompatible with the LDL. 
Whereas the hydrodynamic model has difficulties repro- 
ducing the slope (which may, however, be due to the 
different rapidity windows in the data and the model), 
the LDL gives a completely wrong shape. 

That the RHIC data for semicentral collisions far ex- 
ceed the LDL prediction demonstrates that first collisions 
are insufficient and multiple collisions are required. The 
agreement with hydrodynamics except for a small defi- 
ciency for p t ^1.5 GeV/c indicates that the system is very 
close to local thermal equilibrium. For very peripheral 
collisions we expect that departures from the hydrody- 
namical prediction set in at lower p t -values; it will be 
interesting to see whether experiment confirms this. 

5. Conclusions.- The recent elliptic flow data from 
Au+Au collisions at RHIC show remarkable quantitative 
agreement with the hydrodynamical model, indicating 
a large degree of thermalization in the earliest collision 
stages, well before hadronization. With initial conditions 
tuned to data from central collisions at the SPS and no 
additional adjustment of parameters except for a simple 
scaling of the charged multiplicity to the value measured 
by PHOBOS, the hydrodynamic model reproduces quan- 
titatively the centrality dependence of v-i up to impact 
parameters of about 7fm and its pt-dependence up to 
transverse momenta of about 1.5 GeV/c. Deviations oc- 
cur only in very peripheral collisions and for particles 
with p t > 1.5GeV/c; they may be due to a combina- 
tion of incomplete early thermalization j^] and/or earlier 
freeze-out in these kinematic regions. The low den- 
sity limit LDL roughly reproduces the shape of the cen- 
trality dependence of v-i at RHIC, but slightly underpre- 
dicts the magnitude of the pt-averaged elliptic flow and 
fails badly for the shape of its pt-dependence. It works 
better for the centrality dependence of V2 at the SPS, 
but again cannot describe the observed nearly linear p t - 
dependence of proton elliptic flow. The hydrodynamic 
model gets the shape of all the v-i distributions at the 
SPS right, but seems to overpredict the absolute mag- 
nitude of V2] this last statement is, however, uncertain 
due to the lack of reliable midrapidity data from Pb+Pb 
collisions at the SPS. 

These findings suggest that at RHIC thermalization 
sets in very early (the hydrodynamic simulations point to 
a thermalization time scale of less than 1 fm/c), but that 
it may take longer at the SPS. A better understanding of 
the onset of deviations from hydrodynamic behaviour at 
RHIC, which should be provided by future measurements 
of dN^/dy, V2(pt), and the single-particle p t -spectra as 
functions of the number of participants, will yield crucial 
insights into the kinetic evolution at the earliest colli- 
sion stages. Existing parton |l2],[L3| and hadron M] cas- 
cade calculations reproduce the approximate linear rise 
of v 2 (pt) up to pi^,500MeV/c, but at higher p t the el- 
liptic flow levels off, and in UrQMD its absolute value 
at RHIC, averaged over p t , is underpredicted by about a 
factor 4-5 Jl9j. The parton cascade MPC Jl| builds up 



elliptic flow earlier, but quantitatively does not perform 
very much better. This raises serious questions about the 
adequacy of incoherent scattering among on-shell parti- 
cles to describe the early collision stage and the approach 
to thermalization in ultrarelativistic heavy-ion collisions. 

The accurate agreement of the STAR data || with 
hydrodynamic predictions Q proves that with elliptic 
flow one has found a hadronic signature which is sen- 
sitive to the hot and dense quark-gluon plasma stage be- 
fore hadronization sets in. More detailed measurements 
like those mentioned above (in particular the shape of 
the single-particle spectra) should help to confirm and 
further constrain the picture [ p6| . This will open the 
door to quantitatively characterize the QGP equation of 
state and in particular to distinguish between equations 
of state with and without a phase transition Q . 
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